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Abstract 

The behavior of wake vortices in a convective 
boundary layer is investigated using a validated large 
eddy simulation model. Our results show that the 
vortices are largely deformed due to strong turbu- 
lent eddy motion while a sinusoidal Crow instability 
develops. Vortex rising is found to be caused by 
the updrafts (thermals) during daytime convective 
conditions and increases with increasing nondimen- 
sional turbulence intensity r;. In the downdraft re- 
gion of the convective boundary layer, vortex sinking 
is found to be accelerated proportional to increas- 
ing 77 , with faster speed than that in an ideal line 
vortex pair in an inviscid fluid. Wake vortices are 
also shown to be laterally transported over a signif- 
icant distance due to large turbulent eddy motion. 
On the other hand, the decay rate of the vortices 
in the convective boundary layer that increases with 
increasing 77 , is larger in the updraft region than in 
the downdraft region because of stronger turbulence 
in the updraft region. 
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Nomenclature 

B = aircraft wingspan 

60 = initial vortex separation distance, nB/4 

H = nondimensionai vortex descent distance, 

(z - h 0 ) /b 0 

ho = initial vortex elevation 

L = domain size 

r = radial distance from vortex center 

r c = initial vortex core radius 

S = symmetric part of Vu 

T = nondimensionai time, tVo/bo 

t = time 

t+ = large-eddy turnover time, Zi/w+ 

u,v,w = velocities in x, y, and z directions, 
respectively 

V a = airspeed of generating aircraft 

Vo = initial vortex induction speed, Poo/(27r&o) 

W = aircraft weight 

1/3 

— convective velocity scale, 

[(w , 0 , ) s : surface heat flux] 
x,y,z = axial, lateral, vertical space coordinate 
Zi = depth of convective boundary layer 

zq — surface roughness length 

F = vortex circulation 

To = initial circulation 

Poo = initial circulation at r > r c , 4W/(nBpV a ) 

A = grid size 

€ = turbulence dissipation rate 

?/ = nondimensionai ambient turbulence 

intensity, (ebo) 1 ^ 3 /Vo 
0 = potential temperature 

A 2 = the second negative eigenvalue of the 

symmetric tensor S 2 + O 2 
p = air density 

n = antisymmetric part of Vu 
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I. Introduction 


NASA is developing an automated Aircraft Vortex 
Spacing System 1 (AVOSS) under its Terminal Area Pro- 
ductivity program. This system will determine safe op- 
erating spacing between arriving and departing aircraft 
as based on ambient weather conditions and provide a 
safe reduction in separation of aircraft compared to the 
now-existing flight rules. This system includes prediction 
algorithms 2,3 for transport and decay of aircraft wake 
vortices. To develop this system, research is being fo- 
cused on understanding the interaction of wake vortices 
with the atmosphere. 

A pair of wake vortices usually descends due to their 
mutual induction. The vortex descent rate is reduced 
as the vortices decay with time due to ambient turbu- 
lence and stratification. 3-5 In some cases, however, the 
vortices can rise. It has been shown 0,7 that the vortices 
can either stall or be deflected upwards as they pene- 
trate a strong crosswind shear layer which often occurs 
just above a developing nocturnal stable layer. Vertical 
winds associated with weather phenomena such as vig- 
orous fronts, gravity waves, and thermals may also cause 
the vortices to rise. Rising vortices have been observed 
during daytime convective conditions, 8 possibly due to 
strong thermals developing in the convective boundary 
layer (CBL). Figure 1 shows some observed vertical po- 
sition histories of the vortices during typical convective 
conditions as taken from the Memphis Field Program. 9 
The vortices are observed to remain on the flight level 
(Fig. la), rise (Fig. lb), or even descend faster than an 
ideal line vortex pair in an inviscid fluid, for which the 
sinking rate is given by dH/dT = — 1 (Fig. lc). Numer- 
ical simulation results by Holzapfel et al. 10 indicate that 
in the CBL some segments of the vortex pair that are 
placed on an updraft (thermal) area remain on or even 
rise up to 20 m above their initial positions. As Holzapfel 
et al. pointed out, this situation can be more hazardous 
for reduced spacing operations due to the increased pos- 
sibility of encountering stalled vortices during aircraft 
approach. 

The CBL is typically composed of broad regions of 
gentle downdraft surrounding smaller regions of strong 
updraft (thermal). The fractional areas are about 0.4 for 
the updraft and 0.6 for the downdraft. 11 Both the up- 
draft and downdraft have horizontal scales of roughly 
1.5^, 12 where Z{ is the CBL depth (of the order of 
1 km ). Vertical velocities in thermals can reach 5 m/s or 
more, although weaker updrafts of 1 to 2 m/s are more 
common. 13 Therefore, one can expect a large difference 
in the behavior of wake vortices when they are exposed 
in the CBL, depending on whether they are placed on 
the updraft regions or the downdraft regions. 

In the present study, we investigate the behavior of 
aircraft wake vortices in the CBL using a validated large 
eddy simulation (LES) model. In Sec. II, we describe 
the LES model and the initial conditions. In Sec. Ill, 


(a) MEM-1645 



Figure 1: Time evolution of the observed vertical 
positions of the vortices during daytime convective 
conditions as taken from the Memphis Field Pro- 
gram. 
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we present the results from systematic numerical exper- 
iments in terms of nondimensional turbulence intensity. 
Finally, in Sec. IV, we summarize our results and draw 
some conclusions. 

II. Model and Initial Conditions 

The numerical model used in the present study is 
a three-dimensional, nonlinear, compressible, nonhydro- 
static LES model, called the Terminal Area Simulation 
System 14 (TASS). Numerous successful simulations of 
interaction of wake vortices with the atmosphere 5-7,15-18 
and some within the CBL 19,20 have been conducted with 
this TASS code. In TASS-LES, grid-scale turbulence is 
explicitly computed while the effects of subgrid-scale tur- 
bulence are parameterized by a conventional first-order 
closure model with modifications for stratification and 
flow rotation. 5 21 

The initial vortex system represents a post roll-up, 
wake-vortex velocity field as described in Proctor, 7 and 
consists of a pair of counter-rotating vortices that have 
no initial variation in the axial direction. This vortex 
model has also been adopted for investigating the effects 
of ambient turbulence 5 and ground 10 on vortex decay 
and descent. 

One of the major difficulties in numerical simulations 
of wake vortices is that very fine grid resolution is needed 
to resolve the vortex core whose size is of the order of 
1 m, whereas horizontal domain size of at least four times 
larger than the CBL height (of the order of 1 km) is re- 
quired to simulate large scale turbulence structure in the 
CBL. For current computer resources, it is almost impos- 
sible to simulate large scale structure of turbulence in the 
CBL and much smaller scale structure of aircraft wake 
vortices at the same time. However, the recent study 
by Holzapfel et al. 10 indicates that the main characteris- 
tics of vortex behavior in the CBL can be obtained even 
when greatly reduced domain size (say, about 0.5 km) 
compared to the CBL scale is used. 

In this study, in order to achieve a proper resolution 
of the characteristic scales of both the CBL and the wake 
vortices, we have chosen a domain size of ( L X} L y , L z ) 
= (600m, 400m, 350m) = (155o, 106o, 8.8&o), with a 
grid size of ( Ax , Ay, A z) = (5.0m, 2.5m, 2.5m), where 
y 0 = 40 m, and y and z correspond to the axial, lateral 
and vertical directions of the vortex system, respectively. 
Somewhat larger domain size in the axial direction (155o) 
can allow the development of Crow instability of which 
theoretical maximum wavelength is about 8.65o 22 The 
lateral domain size of 105o is sufficiently large to mini- 
mize boundary influences. The domain size in the ver- 
tical direction (8.85o) is large enough to investigate the 
effects of large turbulent eddies such as thermals on vor- 
tex transport. On the other hand, because of our limited 
computing resources and in order to allow the core to be 
resolved, the chosen initial core size of r c = 5 m is some- 
what larger than the typical value (less than about 5 % 


of the wingspan 15 ) observed behind aircraft. The initial 
vortices are placed on the height of ho = 180m (4.5 bo) 
above the ground. 

Before the vortices are injected, the CBL is driven 
by a constant, spatially uniform heat flux on the bottom 
surface until a typical feature of the evolving CBL is well 
established. Periodic boundary conditions are employed 
in both horizontal directions. Free-slip and no-slip condi- 
tions are imposed at top and bottom boundaries respec- 
tively. The Monin-Obukhov similarity relations are used 
at the surface with surface roughness length, zo = 0.1 m. 
Initial velocities are specified as zero. Initial potential 
temperature is vertically uniform. The turbulence is ini- 
tiated by a small random temperature perturbations su- 
perimposed in five lower layers. 

The vortices are injected at about 3.6 large-eddy 
turnover times (defined as t* = where w* is the 

convective velocity scale), when the CBL is well devel- 
oped although it is still evolving. Figure 2 shows vertical 
velocity fields on the horizontal x — y plane of the initial 
vortex elevation and both vertical and lateral velocity 
fields on the vertical x — z plane along the y— center line, 
with the location of the initial vortex pair indicated. As 
seen in Fig. 2a and b, the right part of the vortices is 
placed on an updraft area with stronger vertical veloc- 
ity, whereas the left part is placed on a downdraft area 
with weaker vertical velocity. From this vertical velocity 
distribution, one can expect that the right part of the 
vortices may stay around the initial vortex elevation or 
rise, striving against the descent by mutual induction of 
the vortices, whereas the descent of the left part of the 
vortices may be accelerated. From Fig. 2c, on the other 
hand, in lower domain layers below the initial vortex ele- 
vation negative lateral velocities are more dominant, and 
thus, the vortices are expected to be transported mostly 
to the negative lateral direction during their descent. 

To evaluate the sensitivity of vortex behavior to am- 
bient turbulence, we generated five values of nondimen- 
sional turbulence strengths of 77=0.20, 0.25, 0.30, 0.35, 
and 0.50 that are in a typical range of the CBL. These 
77 values are typically found during daytime convective 
conditions and obtained by varying the initial circulation 
(Too) to save computing time, rather than by varying ini- 
tial turbulence field. The associated value of turbulence 
energy dissipation rate e is estimated by fitting the theo- 
retical Kolmogorov spectrum in the inertial subrange to 
the simulated spectrum at the horizontal plane of the ini- 
tial vortex elevation, as described in Han et al. 17 Three- 
dimensional wake vortex simulations are carried out with 
each of the five turbulence field, with all other conditions 
identical. 

III. Numerical Results 

Figure 3 shows top and side views of the vortex pairs 
with increasing nondimensional time for three different 
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( a ) w at z= 1 80m 




x (m) 


(c) v along y=0m 



x (m) 


Figure 2: (a) Vertical velocity fields on the horizon- 
tal .r - if plane of the initial vortex elevation and both 
(b) vertical and (c) lateral velocity fields on the ver- 
tical x- z plane along the jy -center line, with the lo- 
cation of the initial vortex pair indicated (thick solid 
line). The contour interval for velocity is 0.4m.?? ! . 


nondimensional turbulence intensities of // — 0.2, 0.3. 
and 0.5. The method of Jeong and Hussain‘S is used for 
the identification of a vortex. They define a vortex in 
terms of the second negative eigenvalue Aj of the sym- 
metric tensor where S and are respectively the 

symmetric and antisymmetric parts of the velocity gra- 
dient tensor Vu. The top (side) view of Fig. 3 is taken 
from x — y (x — z) plane projection of the minimum 
negative A-j value along c (y) direction. A threshold 
value to identify the vortex is arbitrarily determined. 

(a) ETTA = 0.2 


T TOP VIEW SIDE VIEW 



(b) ETA - 0.3 

T TOP VIEW SIDE VIEW 



(c) ETA - 0.5 

T TOP VIEW SIDE VIEW 



Figure 3: Top and side views for the cases of (a) 
r/ - 0.2, (b) r/=0.3. and (c) ;/- 0.5. The unit for axis 
labels is m. 

As shown in the top view of Fig. 3. because of larger 
axial domain size (1 5 ho) used, a sinusoidal Crow instabil- 
ity develops for all three cases, and then, the vortex pair 
connects together as a result of the instability and finally 
disintegrates into a turbulent state. In agreement with 
Han et ales 17 study, the linking of the vortex pair occurs 
earlier with stronger turbulence (that is. larger r/). In 
addition, the vortices are strongly deformed especially 
for stronger turbulence due to large turbulent eddy mo- 
tion. On the other hand, from the side view it is obvious 
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that the right part of the vortex pair rises due to strong 
thermal (updraft) developing in the CBL, whereas the 
descent of the left part vortex pair accelerates due to the 
downdraft. 

Figure 4 depicts the evolution of maximum vertical 
vortex positions that are found in the updraft region. 
It is shown that vortex rise becomes more significant 
with stronger ambient turbulence strength. For the case 
of t/= 0.5, for example, at T= 1.5 the vortex rises up to 
about 2 bo from its initial position and its coherent struc- 
ture is destroyed. For the case of r/= 0.2, on the other 
hand, the vortex appears to stay at its initial height un- 
til about T= 2.5 and then, rises further after its strength 
becomes much weaker due to dissipation by ambient tur- 
bulence. 



Figure 4: Time evolution of maximum vertical vor- 
tex positions for varying nondimensional ambient 
turbulence strengths. 

Figure 5 depicts the evolution of minimum vertical 
vortex positions that are found in the downdraft region. 
It is shown that vortex sinking is accelerated due to the 
downdraft; for all the cases, the sinking rates of the min- 
imum vertical vortex positions are larger than the theo- 
retical sinking rate as given by mutual induction speed 
of an ideal line vortex pair in an in viscid fluid (that is, 
dH/dT = — 1). The sinking rates increase with increas- 
ing ambient turbulence strength; for 77 = 0.6 the sinking 
rate dH/dT ~ —2.3. 

Figure 6 shows the evolution of vertical vortex posi- 
tions averaged along the axial direction and their stan- 
dard deviations. In an average sense, the vortices de- 
scend for all 77 values rather than stay or rise as seen in 
Fig. 6 , although they tend to stall at much later times 
when the vortices almost decay out. This is because the 
downdrafts not only occupy wider area (about 60%) than 
the updrafts, but also the vortex rising due to the up- 
drafts has to counter the effect of sinking due to mutual 
induction. Unlike the behavior of the minimum vertical 



Figure 5: Same as Fig. 4 but for minimum vertical 
vortex positions. 




Figure 6: Same as Fig. 4 but for (a) mean vertical 
vortex positions along the axial direction and (b) 
standard deviations of vertical vortex positions. 
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vortex positions, however, the sinking rates of the aver- 
age vertical vortex positions are smaller than the theo- 
retical sinking rate of dH/dT = — 1 and decrease with 
increasing ambient turbulence strength. On the other 
hand, the standard deviations of vertical vortex positions 
increase with increasing ambient turbulence strength, 
showing about 1.5 bo at T=1.5 for the extremely strong 
turbulence case of rj= 0.5. 

Large turbulent eddies can also transport the vor- 
tices laterally over a significant distance, as indicated in 
the top views of Fig. 3. Figure 7 depicts the evolution 
of maximum and mean lateral vortex displacements from 
the initial positions and their standard deviations for five 
7] values. Although the maximum and mean lateral dis- 
placements and standard deviations appear to increase 
with increasing ambient turbulence strength, the differ- 
ences among varying rj values are much smaller compared 
with those in vertical positions. Figure 7a indicates that 
for the case of r/= 0.2, that is, the most long-lived vortex 
in our experiments, the vortices can be transported up 
to about 2.5 bo at about T = 4 by large turbulent eddy 
motion of the CBL. 

Figure 8 shows the decay of the average circulation 
(averaged over a radius between 5 and 15 m from the 
center) with time for varying ambient turbulence levels 
in both updraft (Fig. 8a) and downdraft (Fig. 8b) regions 
(in each region, the radius-averaged circulations are also 
averaged over 16 planes [80 m] along the axial direction 
where the vortex separation is the widest). As shown 
in Fig. 8, a rapid vortex decay behavior appears to be 
conspicuous due to buoyancy-driven strong turbulence 
typical in the CBL. The rate of vortex decay tends to 
increase with increasing ambient turbulence level, which 
is consistent with field observations 24 and our previous 
numerical simulations. 5 Vortex decay is more rapid in 
the updraft region (Fig. 8a) with stronger turbulence 
strength than in the downdraft (Fig. 8b) region with 
weaker turbulence strength. 

IV. Conclusions 

This study presents the results of numerical simula- 
tions which were carried out for understanding how wake 
vortices behave in the CBL. LES simulations are con- 
ducted by introducing a vortex pair, representing aircraft 
trailing vortices, in a well- developed CBL. Although the 
domain size chosen is relatively smaller compared to large 
CBL turbulence scales, the main characteristics of vor- 
tex behavior in the CBL seem to be well captured. The 
following conclusions can be drawn from the results of 
numerical simulations: 

1) In the CBL, the vortices are largely deformed due 
to strong turbulent eddy motion while a sinusoidal Crow 
instability develops. Consistent with previous studies, 
the linking of the vortex pair resulting from the instabil- 
ity occurs earlier with larger nondimensional turbulence 
intensity. 





T 


Figure 7: Same as Fig. 4 but for (a) maximum and 
(b) mean lateral vortex displacements from the ini- 
tial positions and (c) standard deviations. 


6 






2) Vortex rise is caused by the updrafts (thermals) 
in the CBL, increases with increasing nondirnensional 
turbulence intensity 77 , and can be as large as H ~ 2 at 
T= 1.5 for the extremely strong turbulence case of r/=0.5. 




Figure 8: Same as Fig. 4 but for 5 — 15 m averaged 
circulation in the (a) updraft and (b) downdraft re- 
gions. 

3) In the downdraft region of the CBL, vortex sink- 
ing is accelerated; the sinking rate becomes larger than 
the theoretical value of dH/dT = — 1, increases with in- 
creasing 7 /, and can be as large as dH/dT ~ —2.3 for 

77=0.5. 

4) Wake vortices can also be laterally transported 
over a significant distance due to large turbulent eddy 
motion in the CBL, and their lateral displacements from 
the initial positions increase with increasing ambient tur- 
bulence strength. 

5) Wake vortices in the CBL decay rapidly because of 
strong turbulence intensity typical in the CBL, whereas 
the rate of decay increases with increasing 77 and is larger 
in the updraft region than in the downdraft region. 


These results could serve useful for adding a predic- 
tion of height and lateral variability in the AVOSS vortex 
trajectory forecast. 
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